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Abstract 

In the analysis of the long-time behaviour of two-dimensional incompressible 
viscous fluids, Oseen vortices play a major role as attractors of any homoge- 
neous solution with integrable initial vorticity [5] . As a first step in the study of 
the density- dependent case, the present paper establishes the asymptotic stabil- 
ity of Oseen vortices for slighty inhomogeneous fluids with respect to localized 
regular perturbations. 
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Introduction 

In this paper we consider the motion of a weakly inhomogeneous incompressible 
viscous fluid in the two-dimensional Euclidean space. We can modelize the fluid 
by {p,u), where p = p{t,x) S R"*" is the density fleld and u = u{t,x) S 
the velocity field. The evolution we consider here is governed by the density- 
dependent incompressible Navier-Stokes equations: 

dtp+{u-V)p = 

dtu+{u-S/)u = i(Au-Vp) (1) 
div u — 

where p = p{t, x) € R is the pressure field, which is determined (up to a con- 
stant) by the preservation of incompressibility which yields the elliptic equation: 

div(-Vp) =div(iAit-(u-V)M) . (2) 
p p 
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Alternatively, we can represent the fluid motion using the vorticity field 
10 = curl u G R rather than the velocity. Note that, in the two-dimensional 
context, curl(/i,/2) stands for 81/2 — d^fi- Therefore the evolution equations 
for (p, Lo) become 

r dtp^-{u-v)p = 

\ dtLo + {u-V)uj = div ( i (Vw + V^p) ) ^'^> 

where p is again determined by ([2]) , and u is recovered from uj via the Biot-Savart 
law: 

"(^)-^ / ^p^^(2/)dy (4) 

for X € R^, with (zi,Z2)^ = (—22,-^1)- We also note u = Kbs*^^ where Kbs 
is the Biot-Savart kernel: Kbs{x) = j^jt]^- For simplicity, throughout the 
present paper, the viscosity of the fluid has been rescaled to one. 

We refer to the monograph [TT] for a general presentation of the available 
mathematical results on incompressible Navier-Stokes equations. We also men- 
tion the work of B.Dcsjardins on the global existence of weak solutions |4l[3], 
and, closer to the spirit of the present paper, the work of R.Danchin on well- 
posedness in Besov spaces [1 . Let us emphasize that they both work with the 
velocity formulation ([1]) and do not assume the density p to be bounded away 
from zero. In more physical terms, they allow for regions of (almost complete) 
vacuum, which create technical difficulties. 

In contrast, not only shall we not allow the density to be close to zero but we 
shall only consider weakly inhomogeneous fluids, namely we shall assume that 
the density p is close to a positive constant which, without loss of generality, 
we take equal to one. Remark that if the initial density is contant in space i.e. 
homogeneous, then the density remains equal to this constant for all subsequent 
times. Therefore, in such a case, system ([IJ reduces to the usual incompressible 
Navier-Stokes equations. Moreover, since div(V^p) = 0, the pressure term 
disappears from system ([3]) which thus reduces to 

9t + (w • V) w = A . (5) 

Again, a wealth of information on the Cauchy problem for the homogeneous 
incompressible Navier-Stokes equations can be found in [TT] or [1^. Concerning 
the long-time behaviour of the solutions of the vorticity equation ([5|), the work 
of Th.Gallay and C.E.Wayne has revealed the important role played by a family 
of explicit self-similar solutions, Oseen vortices, given hy p = 1, u — au^ and 
UJ = OLijfi , where 

with 
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and a S R parametrizes this family of solutions. For an Oseen vortex (1, aoj'^), 
\a\ is actually its Reynolds number. If the initial vorticity uq is integrable, it 
is proved in [5| that the corresponding solution of (O converges to aui'^ in L^- 
norm as t — > oo, where a := ^o- Moreover, it was shown in [51 [5] that aw*^ 
is the unique solution of the vorticity equation ([5]) with initial data aSo. 

Even though the homogeneous incompressible Naviers-Stokes equations pro- 
vide a good model in many situations, all real fluids are, at least slightly, in- 
homogeneous and it is therefore important and relevant, both from a practical 
and theorical point of view, to investigate whether the predictions of the ho- 
mogeneous model are meaningful for the density-dependent model, especially 
in the weakly inhomogeneous regime. The goal of this paper is to address this 
question in the particular case of Oseen vortices. These explicit solutions persist 
in the density-dependent case if we assume p = 1, and it is therefore natural to 
ask whether they play there the same role as in the homogeneous case. While 
the general answer to this question is unknown, we treat here one important 
aspect: are Oseen vortices stable solutions for the density-dependent incom- 
pressible Navier-Stokes equations ? In other words, does the theory predicts 
that these self-similar solutions may be observed ? 

Before stating what we mean exactly by stability, let us recall an important 
property of the Navier-Stokes equations: scaling invariance. For any A > 0, 
if {p{t,x),u!{t,x)) is a solution of with corresponding velocity u{t,x) and 
pressure p{t, x), so is 

Dx{p, uj) = (p {\H, Ax), to {\\ \x)) , 

with velocity A w(A^i, Ax) and pressure \^ p{\^t, Ax). Moreover, Oseen vortices 
are self-similar, in the sense that Dx{l,aijj'^) = {l,auj'^), for any a G R and 
any A > 0. To study these solutions, it is therefore more convenient to rewrite 
([3]) in self-similar variables. Following [7] , we set 

(r,0:=(lni,-^) . (6) 

Motivated by scaling invariance, we will work with new quantities {r,w,v,lVj 
related to the former by 



Then the corresponding evolution equations for (r, w) are 

w + ((w - i ■ V) u; - w = div ( i (Vt^ -I- V^H) ) 



(7) 



(8) 



where again v is obtained from w by the Biot-Savart law and VII by solving 

div ( i Vn ) = div ( i A t) - (w • V) v) . (9) 
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By constuction, Oseen vortices correspond to stationary solutions of ^ , namely 
(1, a G), where a is any real number. We fix initial data for the original equa- 
tions a,t t — 1 rather than at i = 0, hence at r = for the new equations. 

In order to state our main result, we now write down the evolution equations 
for a perturbation of an Oseen vortex. Let a € R. We work with new quantities 
{b, w) related to the former by 6 = r ~ ^ ^^'^ w = aG — w. The evolution 
equations for (6, w) are 



drb+{{v^^O-V)b = 
drW - (C- aA)w + {v ■\7)w = div ( 5 (Vw + V-^H) [ 



(10) 



where 

cf = A/ + i^v/ + /, 

A/ = t;G.V/+(XBs*/)-VG 



(11) 

Here v is obtained from w by the Biot-Savart law, w and v are recovered by 

w — aG + w , V — av^ + v , (12) 

and Vn is obtained by solving 

div((l + 5)Vn) =div((l + &)Aw- (w- V)u) . (13) 

Let us point out that, thanks to the linearity of the Biot-Savart law, we do have 
V = Kbs * w and v — Kbs * w. 

Note that C is the usual Fokker-Planck type operator which generates the 
evolution corresponding to the heat equation in self-similar variables. On the 
other hand, A is the non-local first-order operator resulting from the lineariza- 
tion around w — G oi the transport term of ([8]). More precisely, if = Kbs*w 
and {v, w) satisfies (fT2l) . then v ■ Vw = aAw + v ■ Vw, since v'~^ _L VG. 

Before stating our result of stability, we finally introduce the function spaces 
and norms we will encounter throughout the present paper. For 1 < p < oo, we 
denote by L^'(R^) the usual L^'-space and by | / |p the usual i^-norm of a func- 
tion /. Similarly, for s G R, we denote by iJ*(R^) the usual Sobolev space and 
by I / the corresponding norm. Sobolev norms will be convenient to specify 
the regularity of our perturbations, but we also need weighted norms to describe 
the localization properties. Indeed, even in the homogeneous case, we know that 
it is impossible to get a convergence rate in time if we do not assume that the 
perturbations decay sufficiently fast at infinity in space (see [8]). Instead of 
using polynomial weights as in [S], we choose here the Gaussian weight G~2^ 
which is naturally related to the Oseen vortices and has several technical advan- 
tages. For instance, on the Hilbert space L^(R^) := {/ | G~2 / g ^^(R^)}, the 
linear operators C and A become respectively symmetric and skew-symmetric. 
Consequently, for any 1 < p < -t-oo, we shall use the weighted L^-space defined 
as follows: 

LPJK') :^ { / I G-i / e LP{R') } (14) 

with the corresponding norms | / \^^p :— \ G^i f \p. 

We are now able to state the main result of this paper: 
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Theorem 1 Let a G R, < s < 1 and < 7 < ^. There exist Eq > o.nd 
K > such that for any < e < £0 */ 

1. bo e iJ^+2(R2) is such that G^Hq belongs to L2(R2) n L°°(R2) with 
I ^0 \w,2 < £, I &0 \w,oo < ff^c? G^2 V60 to L''(R^) /or some q > max(4, |) 

^. W() E H'^CR?) is such that G^^wq belongs to L^(R^) with \ 'Wo\w,2 < e, 
and Jj^2 Wo = 

i/ien t/iere exists a unique solution (6, 0/ U0\) with initial data {bo,WQ) such 
that 

1. 6g L^,(R+;H^+2(R2)), 

2. G-ifee L°°(R+;i2(R,2)p2.°°(R2)), G"* V5 e Lj^^(R+; ^'(R^)) 
5. w G L^,(R+;i?^(R2)) n Ll^{B.+ -H^+\R^)), 

4. G-^w e L°°(R+; L2(R2))nL2(R+; L2(R2)), G-^ G L2(R+; L2(R2)), 
G-^\^\weL^{R+;L\R^)). 

Moreover this solution satisfies \ w{t) |u,^2 < K ee~^'^ , for any r > 0. 
Let us make some comments on this result. 

1. We treat only small regular and localized perturbations. Note also that 
one can not allow for perturbations with non-zero mean vorticity. Indeed, 
since Jj^a w is preserved under the evolution, in order to have convergence 
in L^(R2), we must have a = wq. Consequently, henceforth, we will 
scarcely repeat but always assume wo = 0. 

2. In contrast, one important aspect of our result is that we make no small- 
ness assumption on the size of the Oseen vortices, so we do treat high 
Reynolds numbers |a|. Yet note again that the present paper is far from 
telling us whether Oseen vortices attract all weakly inhomogeneous fluids. 

3. Now, to make our result more concrete, let us give some immediate conse- 
quences in the original variables. Under the asumptions of Theorem [l] if 
{p, Lu) is the solution of ^ defined by ([7]) with w = aG + w and r = 
then the vorticity uj(t, x) satisfies 

ti\e^ iujit)-au;^it))\2 < ^ , t>l. 
Moreover this implies 

ti-^|w(t)-a^«(t)|p < ^, ti--.\uit)-au''{t)\, < ^ 

1 _ 1 1 — 1 

for any I < p < 2, 2 < q < +00, t > 1. Note that t p and t^ 1 are 
the rates of convergence corresponding to Oseen vortices. 
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4. Concerning the density, we can not expect the perturbation to decay in 
time in the whole space, since we deal with an incompressible fluid. Indeed, 
since the density satisfies a transport equation by a divergence- free velocity 
field, formally it remains constant in law. Thus any L^'-norm of a function 
of p remains constant in time. 

To establish Theorem [U following , we build a sequence of solutions of a 
linearization of (jlOp . In order to show the convergence of this sequence, we will 
use local-in-time estimates in Sobolev norms for solutions of this linearization. 
To be more precise, in a first step, we establish global estimates for solutions of 
the linearized system, which control the density in weighted U' spaces and the 
vorticity in . In a second step, we use the previous estimates to prove first 
global estimates of the density in H'^~^ (with a loss of regularity), then local 
estimates of the vorticity in H'^^^ , and finally local estimates of the density in 
ffs+2 j^i^Qj^ these estimates to establish the existence and uniqueness 

parts of the theorem. Finally, we prove the asymptotic part of the theorem. 

The plan of the paper does not follow the former steps, but we will encounter 
the structure of the proof in the various sections. In a preliminary section, we 
collect estimates on the Biot-Savart kernel, thus on the velocity in terms of the 
vorticty, and on the pressure. In the second section, we gather estimates on the 
linearized density equation. In the third one, we establish the estimates on the 
linearized vorticity equation. At last, we state and prove our main results. 

In what follows, the original (and physical) time will never appear again, so 
for notational convenience henceforth we allowed us to use the letter t to denote 
the rescaled time r. 



1 Preliminaries 

If / is intcgrablc over , we define its Fourier transforms to be the function / 
defined for any rj G by 

fiv)= f f{C)e"^<dC . 

Concerning function spaces, we will also need the following convention. For 
any cr S R, we denote by 7j'^(R^) the usual homogeneous Sobolev space on R^ 
equiped with | / 1^,, := | /'^/ I2, where / := (— A)2 . 

1.1 Biot-Savart kernel 

The goal of this subsection is to enable us to estimate the velocity in terms 
of the vorticity. We collect here some estimates on v in terms of w when v is 
obtained from w by the Biot-Savart law, thus on the velocity in terms of the 
vorticity. Recall that, in this case, for almost every ^ e R^, 
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where {xi,X2)'^ = (—X2,xi), that is v = Kbs * w, where Kbs is the Biot- 
Savart kernel Kbs{x) = "i^f^^- Also note that in terms of Fourier transform, 
(fT5|l becomes 

^ ■ (16) 

Most of these estimates are already known, but for the sake of completeness we 
try to give proofs rather than references. 

The following proposition gathers estimates in L^'-spaces. 

Proposition 1 

1. Let 1 < p <2 < q be such that 1 + - = ^ + -. 

^ ^ p 

There exists C > such that, if w belongs to LP(R^), then defines a 
V in L'^(R^) and 

\v\,<C\w\p. (17) 

2. Letl<p<2<q< +oo and < 9 < I be such that ^ + = i 
There exists C > such that, if w belongs to LP(R^) n L'^(R^), then fZ5)J 
defines a v in L°°(R'^) and 

\v\^<C\w\l\w\l-'^ . (18) 

3. Let p > 1. There exists C > such that, if w belongs to LP(R^) and v is 
defined from w by US]) , then Vw belongs to L^'(R^) and 

\Vv\p<C\w\p . (19) 
In addition, in all cases, we have divv = and curlu — w. 

Proof. Part 1 follows from a Young-like inequality called Hardy-Littlewood- 
Sobolev inequality (see for instance Theorem V.l in [12] )■ Indeed, Kbs is 
weakly L^ but not square intcgrable. 

Part 2 is trivial when w = 0. If not, we remark that from Holder's inequali- 
ties, we obtain 

< C \ w\g R'-l +C \ w\p ^ , 

Rp 

for almost every ^ S R^ and any R > 0. Aiming to optimize this inequality, we 
choose R = O-nrr)'^ with /3 = 4^ = and derive d]). 

Part 3 holds since derivating ([TS]) yields that Vw is obtained from w by a 
singular integral kernel of Calderon-Zygmund type (see Theorem II. 3 in [H]). 



The following proposition gathers estimates in Sobolev spaces. 



7 



Proposition 2 

1. Let s e R. There exists C > such that, if w belongs to _ff''~^(R^), then, 
if V is defined by v belongs to _ff*(R^) and 

\v\fj.<C\w\fj.., . (20) 

2. LetO< s < 1. 

There exists C > such that, if (1 + |^|) w belongs to L^(R^), then, if v 
is defined by 173)) . w belongs to H'^{'R'^) and 

\v\H.<C\il + \^\)w\2. (21) 



Proof. Part 1 is a direct consequence of the Fourier formulation (fT6|) . 
Part 2 is thus reduced to estimate | it; . If < s < 1, we note that 

The second term of the right member is dominated by | w ||. Choosing p such 
that P > I and applying first Holder's inequalities then Sobolev's embeddings, 
we obtain for the first term 

Finally, gathering every piece yields 

\w\%^, < C\w\l,+C\w\l < C\{l + \^\)w\l . 
This concludes the proof. ■ 



1.2 Pressure estimates 

Keeping in mind equation p3p . we gather some estimates for a solution 11 of 
the following equation: 

div((l + 6)Vn) =divi^ . (22) 
We begin with estimates in L^-spaces. 
Proposition 3 
1. Let p> 1. 

There exists C > and k > such that if F belongs to LP(R^) and b 
to L°°(R^) with k|5|oo < 1, then \22\) has a unique solution 11 (up to a 
constant) such that VII belongs to L^(R^), and 

|Vn|,<^-^|^|,. (23) 
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2. Let 1 < p,q,r be such that :^ = ^ + ^ • 

There exists C > and n > such that if F belongs to L'^(R^) and, for 
i = l,2, hi belongs to L°°(R2) n LP(R2) with 

I |oo ^ 1 o,f^d Tli solves 
div ((1 + 6i)Vn,;) = divi^ , 
then V (112 — Hi) belongs to L''(R^) and 

|v(n2-ni)|, <-— ^i_|62-6i|p|i^|, . (24) 



Proof. 

1. We want to obtain VII, the solution of , in terms of F as a perturbation 
of Leray projectors. Let P be the Leray projector, that is the projector on 
divergence-free vector fields along gradients, and let Q = I — P. Remark 
that Q F gives the solution VII of ((22|) when 6 = 0. Now we can rewrite 
([22]) as 

divVn = div(F-5Vn) 

then Vn = Q (i^ - 6 VH) thus (I + Q 6) VH = Q i^. Since Q is continuous 
on LP, there exists k > such that 

IQbflp < k|6UI/Ip 
for / e LP(R2). Thus, if K I 6 loo < 1, I + Q 6 is invertible on LP, and 

Vn= (I + Q6)-iQF (25) 
gives the unique solution, with the expected bound. 

2. Reminding (|25p we write 

VHi = (i + Q&2)-Mi + Q&2)(i + Q&i)-'QJ^ 
vn2 = (i + Qfe2)"Mi + Qfei)(i + Q&i)"'Q^ 

then subtracting and factorizing yields 

V(n2 - Hi) = (I + Q 62)"' Q {bi - 62) (I + Q &i)-i Q F . 

Now the continuity of the operator Q on L^ reduces (|24|) to an estimate on 
I (61 — 62) (I + Q ^i)~^ Q F \r- At last applying first Holder's inequalities 
then the continuity on concludes the proof. 



In order to estimate solutions of (|22|) in Sobolev spaces, we first state some 
useful commutator estimates of Kato- Ponce type (see [^). Let us recall that 
/ = (-A)i. 
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Lemma 1 Let < s < 1 and a > 1. 

1. There exists C > such that if I""/ belongs to L^(R^) and g to i7'^(R^), 
then P{fg) ~ f I^g belongs to L^(R^) and 

\lHf9)-fI'9\2<C\rf\2\g\H^. (26) 

2. There exists C > such that if I'^f belongs to H'^{R?) and g to L'^(R'^), 
then P{fg) — f I^g belongs to i^(R^) and 

\I\f9)-fI'9\2<C\Pf\H^\g\2. (27) 



Proof. Let us first note that there exists C > such that 

l^ll < C\h\H^ , 

for any h in iJ°'(R^). This comes applying Holder's inequalities to 

(1 + byr)^ 

Therefore in order to prove the lemma it is sufficient to establish 

\I%f9)-fI'9\2 < C\Pf\2\9\i, (28) 
\I'if9)~fI'9\2 < C\I^f\,\g\2. (29) 

Now set h = r{fg) - f Pg. We have 

'^('?) = Wxi I (1^1' - 1^ - ^1') ^(^ - dC 

for almost every rj E R^. Thanks to the foUowing basic fact: 

M'-Wn < h-ryT, 0<5<1 (30) 
for rj,ri' E R^, we obtain 



I lCn/(C)l \9iv-0\dC 



At last depending on how we apply Young's inequalities we obtain either 
or (Hi). 



We now state the anounced estimates in Sobolev norms. 
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Proposition 4 Let < s < 1 and a > 1. 

There exists C > and k> such that if F belongs to i/-'(R^), b to L°°(R^) 
with K\b\oo < 1 and rb belongs to H^iR'^), then, if U solves (MW, VH 
belongs to L^(R^) and 

U^Vnb< ^ (|rf|2 + - \j;^\I'b\H^\F\,) . (31) 

Proof. Applying to (|22p and commuting 6 and yields 

div ((1 + b)\/P n) = div ([6, /"] Vn) + div (/" F) . 
Applying then (|23p to this equation we obtain 

|/^vn|2 < - — ^1— (|rf^|2 + |[6,r]vn|2) . 

1 - K I &|oo 

Now using first ([ST]) then applying ((23)) once again yields 

\[b,r]wu\2 < c\rb\H^\wn\2 < - — ^T^irM^H^h • 

1 - K I loo 

Gathering everything leads to ([5T|) . ■ 



2 Density equation 

In this section, we gather information on the following linearization of the den- 
sity equation: 

ft 6 + ( (i. - i e) ■ V ) & - (32) 

where z7 is a divergence-free vector-field, a € R and v = av^ + 1^. By lineariza- 
tion, we mean that we do not assume that v is obtained from a solution w of 
the vorticty equation in (jlOp . which involves b. 

We begin with estimates in L^-spaces, weighted or not. We recall that 
L^(R^) for 1 < p < 00 is the weighted space defined in (fT4|) . The next propo- 
sition is the density part of the annouced first step of the proof. 

Proposition 5 Let T > 0. 

Assume that v is a divergence-free vector field belonging to L^(0, T ; L°°(R^)). 
Then b, the solution of \3S]) with initial data bo, satisfies 

1. for any \ < p < +00, provided 69 G L^i^'^), 

\b{t)\p <\bo\pe-i , forO<t<T- (33) 
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2. for any I < p < +00, provided G 2 5o G L^(R^), 

< |&oU,pe-pe^/o forO<t<T. (34) 

Proof. 

1. For 1 < p < +00, multiplying ([32|) by sgn(6) where sgn is the usual 
sign function, and integrating by parts yield 

since diviy ~ and div^ = 2. Integrating gives (|33p in this case. The case 
p = +00 follows letting p go to infinity. 

2. For 1 < p < +00, starting as in the former case leads to 

i I ' 1-^ — L ( ('^-^ ) 1^1' - L ( i ^•('^-^ ) 1^1'' 

since VG^^ = | ^. Now, since ^ _L f'^(C), we have 
hence 

||6|S,,p<(-i + |l^^lL) 

Again integrating achieves the proof for finite p and the case p = +00 
follows letting p go to infinity. 



The next proposition corresponds to the density part of the announced sec- 
ond step: local-in-timc estimates in Sobolev norms. In order to prove some part 
of it we will need the following commutator lemma. 

Lemma 2 Let s > 1 and a > 1. 

There exists C > such that if I"^ f belongs to L^(R^) and g to H'^{'R?), then 
I'ifg) - f I'g belongs to L^(R^) and 

I I'ifg) -fI'9\2<C\ Ff I2 I 5 lif^ + C I V/ I P-'g b (35) 

Proof. The proof is essentially the same as the one of Lemma 1 except that 
here ([50]) is replaced by 

Iicr-ICTI < C|C-C'I(IC-CT"' + ICT"'), forCC'GR'. (36) 
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Note that we could have obtained, as in Lemma 1, various estimates depending 
on the way we apply Young's inequalities. ■ 



Proposition 6 LetT >0. 

1. Let0<s<2 and < e < s. 

Assume that v is a divergence-free vector field with V? e i^(0, T; _ff^(R^)). 
Then there exists Ct > independent of v such that, for any initial data 
bo e i7"(R2), any solution b G L°°{0,T; H'-'(R'^)) of (3^ satisfies 

I bit) \hs-. < Ct I 6o \h^ exp ((Ct / | Vi/(t) j^i drf) , forO<t<T. 

Jo 

(37) 

2. Let s>2. 

Assume thatv is a divergence-free vector field whith "S/v G i^(0, T; H'^~^{T{?)). 
Then \32\) has a unique solution b G L°°{0, T; if*(R^)), for any initial data 
bo E H^{Tl^). Moreover there exists C > independent of bo and v such 
that b satisfies 

I U= < 1 5o U= e^* exp (C / \Vv{t)\h^-i dr) , forO<t<T. 

Jo 

(38) 



Proof. 

1. See Theorem 0.1 in [2]. 

2. We can compute §] ' / ~ /^^^ div/, for any vector field /. Thus 
applying /" to equation p2p and commuting yield 

9t/^^>+((i'-io•v)r6 = ^r6-[/^H• V6. 



Then multiplying by Pb and integrating lead to 

1 d , ,o s 

2 dt 



since div ly = 0. Now use Cauchy-Schwartz' inequality and apply Lemma 2 
(with (T = s — 1) to get 

At last combine the former with -^^Ifcll + ^l^li — Oto obtain 
which yields (j38p by a mere integration. 
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The last estimate we state for the linearized transport equation (l32|) is in- 
tended to be used for the proofs of the convergence of our iterative scheme and 
of the uniqueness of our solutions. Indeed we estimate the difference of two 
sohitions of equations of type (|32p . 

Proposition 7 LetT > 0. 

Assume that, for i = 1,2, is a divergence-free vector field belonging to 
L2(0,T; Vl^i'°°(R2)). //, fori = 1,2, 6, ts a solution of 

where Ui — av^ + T'i, with initial data bo, then bi and 62 satisfy 

1. provided that G^^Vfeo belongs to _LP(R^), for some 1 < p < +00, 

|V5.(t)k, < |V6okpe-*(t-^)e*/o l^-(^) Jo IV''.(«)U<^« (39) 
/or i 1,2 andO <t< T; 

2. provided that bo is such that G^^bo belongs to Lp(R^) and G~^V6o to 
L^(R^), for some 1 < p < q < +00, 

|(62-6i)WU,P <e^^« sup |V6i(s)U,, / \{i;^-i;,){s)\rds 

0<s<t Jo 

(40) 

for < t < T , where r is such that ^ = ^ + 
Proof. 

1. Derive the equation for bi to get for j ~ 1, 2, 

dt djbx + d,bi = -djvi • V61 + i djbi . 

From this, following the proof of (p4|) . we obtain for j = 1,2, 

|| dMt) < (-1+1+1 1 ^1 W I-) I dMt) \l.v+p I V'^iW U I v&iW 1^ 

Now combining the incquahties for j = 1,2 and integrating lead to p9p . 

2. Observe that 62 — satisfies 

dt (62 - 61) + ( (i^2 - ^ • V) (62 - foi) = -(i?2 - ^1) ■ V61 . 
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Now following again the proof of yields 




where (5 = 62 — fci- It is now straightforward to derive (PD|) . 



3 Vorticity equation 

As announced, we now study a linearization of the vorticty equation: 

dtw - (C ~ a A) w + ■ V) w ^ diy {b (Vw + V^H)) (41) 

where C and A are as in (jlip , 6 is a real function, v is a divergence-free vector 
field, a e R, 

V — Kbs*w , V — av^ + v , 
w — aG + w , V = aw'^ + z?, 

and Vn is obtained by solving 

div ((1 + h) vn) = div ((1 + 6) Au - (i/ • V) v) . (42) 
Remind that we always assume /^^a = 0. 

3.1 Global estimate 

In this subection we establish a global-in-time estimate in weighted L^-spaces. 

Proposition 8 Let a £ R, Kq > 0. There exist £0 > and C > such that if 
b is a real function and v a divergence-free vector field such that 

1. for < t < T , for any 1 < p < +00, 2 < q < +00, 

2. forO<t< T, 

3. and 

I ^0 \wA < £0 : \bo \w,oo < £0 
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then any solution w G L'^(0, T; L^(R^)) of [^\ ), with initial data wo G L^(R^ 
satisfies, for any < t < T, 



w 



it) 11,2 + C [\\w \l, + I \l, + I \^\w (43) 
Jo 

< 2|7?oLl2 +C|a||6oU,4. (44) 

Note that the assumptions on b corresponds to the first proposition of the 
previous section. Note also that once a and Kq are fixed, since is embedded 
in any L^, 1 < p < 2, and is embedded in any L'', 2 < q < oo, inequahties 
and ([33]) enable us to make | w |^ as small as we want provided we take 
Wo and bo small enough. At last note that the proposition enables us to bound 
lo I lifi ' which wiU be used in ((57)) . 

Proof. Our strategy is to multiply (|4ip by w and integrate to bound 
'm\'^\w 2- what follows, we examine each term arising once multiplyed by 
w and integrated on R^. 

• Let us emphasize first that (j41]) preserves Jj^a w. Hence j-^^ w = 0. 

Let L :— G^^ {—C)G^. A direct calculation shows that L = — A+ \ 
is a harmonic oscillator with spectrum {0, i, 1, f , • • ■}• Moreover is a 
simple eigenvalue with eigenvector G^ . In particular, if / belongs to the 
domain of L with Gi f = 0, then f^f> ^ I / ll- 

Coming back to C, we obtain: if G^^w belongs to the domain of L with 
/r2 w — 0, then, for any < 7 < ^, 



G-^w Cw < --{l--f)\w\l2 + -f I G-^wCw 



1 

/R2 2 J^2 

thus integrating by part from the formula for L 

/ G-'wCw < -Ul~2j)\w\l,~^{\ViG-iw) 

and expanding 
S^,G-^wCw < -i (1-27) Iwp 



w,2 

lii,7;|2 I 9 r r:-^x7-,r,. L 



hence 

G-^ w Cw < [1 ~ 2j) \ w\l ^ - \ Vw \:„ ^ + - \ ^ w 



R2 



Recalling that Kw — v'^ ■ Vw + v ■ VG, we obtain 



R2 



(45) 

7G, we obtain 
G-^wAw = 0. (46) 
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Indeed , from u'^(C) -L (, and VG ^ = — |G we derive 

And, on the other hand, using the identity r]-^ = — -rj and the expHcit 
formula (fT5|) for the Biot-Savart law, we derive 

= Jr2xr=^(0 w{ri)dr)d£, 

= - J^,G-^wv-VG . 

Thus G-iww-VG = 0. 

Using Holder's inequalities, we also obtain 

\ j^^G-'wv-Vw\ < 6\^^\l\w\l,2 + ^\^iS\l^^ . (47) 

Integrating by part, we obtain 

/ G-^ w div (bVw) ^ ~ [ G-H\Vwf~l[ G-^bw^-Vw 

JR2 JR2 2 7r2 

then, using Holder's inequalities and the fact that | b{t) |oo < | 6o |oo, 
I J^^ G-i w div ibVw)\ <l\boU\Vw + i I &o loo I 1^1 11,2 ■ (48) 

In the same way, since | b{t) I2 < | &o be"^ , we have 

|^.G-i^;div(6VG)| = \y^.^b( -Vw + l J^,bw\^\^\ 

< Gl^obe-* (|V?UU,2 + |wU,2) 

thus 

1/ G-iz? div(6VG)| < G|6o|2(e-* + |ViC|S,.2 + |55|^,2) • (49) 

JR2 

Finally, integrating by part, using Holder's inequalities and applying in- 
equality we obtain, for bo small enough in L°°, 

\[ G^i w div(6V^n)| = |i [g-HwC-V^U+[ G^HVw-V^ 

JR2 2 J J^2 

< 1 ^J-^{\\i\w\^,2 + \^w\^,2) 

1 - K I 60 loo 

X (l^oU.ooe^" |(H-5)AiJ|2 

+ I bo 1^,4 e^" e-i\a{l + b) Av^ - (i^ • V) w 
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Now, on one hand, since | b{t) |oo < | ^^o |oo, estimate (PU]) yields 

\{l + b)Av\2<C{l + \boU)\Ww\2 . 
On the other hand, similarly, we have 

1(1 + 6) At;« |4 <C(l + |5oU) . 

At last. Holder's inequalities, estimate p9|) and Sobolev's embeddings 
yields 

\{iyV)vU<C{\a\ + \I^\s){\a\ + \w\s)<C{\a\ + \I)\s){\a\ + \w\Hi) . 
Taking this into account yields when k | feo |oo < 5 

I / G~^w X div(6V-Ln)| 

< e-i C|a|e^'M&oU,4(l + |a| + |J;|8)' 

+ 1^^11,2 Ce^°|6oU,4(l + |a| + |^7|8)2 

+ |Vi^|^^2 Ce^«(|6oU,oo + |fooU,4(l + |a| + |i?|8)) 
+ lk|5^l^,2 C^e^«(l&oU,oo + |6oU,4(l + |a|)) . (50) 

It only remains to us to gather everything after setting 7 = i in (|^ and 
integrate in time in order to obtain, when k | 60 |oo ^ 



w{t)\l, X (l-12/„*|i;|^) 










X 


(1- 






X 


(1- 


- 


+ cj;^\\^\w\i^. 


X 


(1- 


- e^'> 1 




X 




\w,4 6 



which yields the proposition since /q | | ^ — 



24- 



3.2 Sobolev estimate 

In this subsection we prove a local-in-time estimate in Sobolev norms for solu- 
tions of equation (|^T|l . Remind that / = (— A)2 . 

Proposition 9 

Let < s < 1, 1 + s < s' < 2, 1 < s" < 2 - s and a G R. 
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There exists Eq > and, for K > 0, there exists C > such that if b is a real 
function and v a divergence-free vector field such that 

sup[o,T] I ^loo < |&o|oo<eo , supjo y] l&l^,' < K , 



/o'lV^ili < K , /; < K , 

then any solution w € L°°(0, T; ij*(R^)) of |77P , wif/i initial data wq e ij^(R^), 
satisfies for any < t < T, 

\rw{t)\l + C f\I'Vw\l < C e""' {\Pwo\l+K) . (51) 
Jo 

Note that ([551) ^'^'^ ([57)) could provide us the validity of the asumptions on 
b, and (j43p both the validity of the asumptions on w and, thanks to estimates 
([T5|) , (I^Dl) and d^T]) , the validity of estimates on V when v = v. Conversely (|5ip 
could be used in ()38p . again when v = v. 

Proof. We choose a such that 1<(t<1 + s,s + (7<s' and a < s" . The only 
role of a is to make clearer our use of commutator estimates ([26]) and ([27]) . 

Our strategy is to apply to ([¥T|) . then multiply by /^u; and estimate each 
term arising, in order to bound ^| Pw |^2- 

• First we compute the commutator — ^I" and obtain 

/ rw rcw = - f \vrw\^ + f ir^p (52) 

since integrating by parts yields /j^^ f Cf ^ - J^^ I V/p + i /j^^ 

• On one hand, we have 

rw r{{v^-v)w)= f rw iv^-v)rw+f rw i[r,v^]-v)w 

R2 JR2 Jr2 

with, integrating by parts, 

/ Pw {v^-V)rw^lf (v<= • V) = 

JR2 2 Jr2 

and, using Holder's inequalities and inequality ([27]), 

1/ r^; i[r,v'']-v)w\ < c\rw\2\rv''\H'\yw\2 . 

JR2 

On the other hand, we have 

/ rw /"((?• v)G) = / rw {v-\7)rG+ [ rw {[r,'U\-V)G 

■/r2 JR2 Jr2 



19 



with, using Holder's inequalities, 

1/ rw iv-v)rG\<c\rw\2\vW^rG\2 

and, using both Holder's inequalities and inequality (1261) . 

1/ rw {[I^^-^)G\<C\rw\2\I'v\2\\7G\H^ . 

Using estimate (PT|) and estimate ([T5)) combined with Sobolev's embed- 
dings, we derive 

I / rw PAw\ < C{\w\l^^ + \Vw\l) . (53) 

• In the same way, we have 

/ Pw (?-V)rw==0 
and, applying Holder's inequalities and inequality (|26p . 

I / PVw-i[r,^w)\ < G\rVw\2\P^\2\w\H^ 

JR2 

thus, since < cr < 1 + s, 

1/ Pw r{{i^-w)w)\ < G\ri:^\2\rww\l + G\rv\2\w\l . (54) 

JR2 

• Integrating by parts yields 

/ rw rdiv(5Vw) - / b\r\/w\^~ [ r\7w [r,b]\7w. 

Applying (|27p . we derive 

I / rw rdiv{ bww)\ < (I bo \oo + e)\ rww \l + -\b | ww \l 

(55) 

where e > is intended to be chosen small enough. 

• Similarly, 

1/ Pw rdiv(6VG)| < (|6o|oo+£)|/^V^;|2+C(|6o|oo + -|&llf=+.) • 
Jr2 £ 

(56) 
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• First, integrating by parts yields 

/ rw /Miv(w-Ln) = - / rvw-brv^n-f rvw-{[rMv^ii) 

with, using Holder's inequalities, 

I / FVw ■ brv^ii\ < |rvw|2|6o|oo|/'vn|2 

and, using Holder's inequalities and inequality (I27p . 

I / rVw ■ ([r , 6]V^H)| < c I rvw u \ b \hs+^ I vh I2 . 

Besides, on one hand, estimate ([23)) applied to equation (|42l) and Holder's 
inequalities imply 

IVHI2 < Ci\a\ + \a\^ + \Av\2 + \a\\\/v\2 + \uU'^v\2 + \a\\uU 

thus with estimates on Biot-Savart kernel 

IVHI2 < C(|a| + |a||I7|oo + |a|2 + (|a| + |I7|oo)|S^|2 + |Vw|2) . 

On the other hand, estimate ([5T|) applied to and Holder's inequalities 
imply 

|rvH|2 < c{\b\Hs+.\vu\2 + \r{{i + b)Av)\2 + \r{iyV)v\2) 

with, commuting P and b thanks to ([27]) . after some calculation, 
\riil+b)Av)\2 < C\b\Hs+.i\a\ + \Vw\2)+Cil+\bo\oo)i\a\ + \rVw\2) 
and, in the same way, commuting and i', 

\imiyV)v)\2 < C{\a\ + \I'^\H'){\a\ + \il^\2)+C{\a\ + \l)U{\a\ + \rw\2). 

Therefore, gathering these inequalities, 

\J^2l"w X /''div(6V^H) I 

< C\rw\l\bo\^{\a\' + \l^\'^) 
+ C\rVw\U\bo\oo+e) 
+ C|6|^.+„(i + |6o|oo) 

x{\a\' + \a\^\V\l + \a\^ + {\a\^ + \i;\l)\w\l + \yw\l) 

+ C|5o|oo(|ap(l + |6o|L + l«P + |i^Poo) 

+ + I ) (lap + I 55 li) + I 6 11,.+. + I li) ) . 

(57) 
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Putting all these points together and integrating yields 



\rw 


1 

2 








(l-|6o|cc^(l + |a|)-sup[o_,]|/«?|2) 


< 1 1'wo 


\l 




+ cf. 


rw\l 


(l + |6o|oo/o|i?|L + |a|(l + |6o|oo/ol^^lLl«||6o|oo)) 


+ 




(l^^oloolap 






+ + sup[o \w\l){l + \ bo loo) sup[o_t] 1 b ) 


+ ^/o 


Vw\l 


(|a| + (l + |6o|oo) sup[o_t] |6|^,+„) 


+ C/o 




( |a| + lap 1 6o loo + sup[o,t] 1 1'^ b 






+ |a|2(l + |6o|oo) sup[o,t] |6|^,+„) 


+ C^/o 




1^0 loo (lap +sup[o,t] 1 w|i) 


+ C t 




|a| (1 + lap) (1 + 1 6o loo) (1 + sup[o,i] | b |2^.+„) . 



(58) 



A Gronwall-type argument achieves the proof. ■ 



3.3 Estimate for convergence 

We now establish an estimate on the difference of two solutions of equations of 
type (|4T|) , intending to prove convergence of our iterative scheme and uniqueness 
of solutions of (fTU|) . 
For i = 1,2, consider 

^tw^-{C~aA)w^ + {'i:^^■W)wi^ div {bi (Vw, + V^n,)) (59) 
where C and A are as in (jlip . hi, fli are real functions, a G R, 
Vi = Kbs *m , h'i = Kbs * , 

Wi — aG + Wi , rii — aG + fli , 
and Vlli is obtained by solving 

div ((1 + b^) Vn,) = div ((1 + bi) Aw, - (z^, • V) V,) . (60) 

Note that we choose to write — Kbs * ^ to stress the symmetry of the 
hypotheses on fl and w. 

For concision's sake, we note (5/ = /2 — /i, for any functions /i, /2. 

Proposition 10 

Let a eH, K > 0, (T > 2, < T] < s < 1 and max(|,4) < p < +oo. 

There exists £o > and, for K',T > 0, there exists C > such that ifiHi, W2 

satisfy / f?3)) with 
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1- I bo \w,4 < £0 , \bo \w,oa < £0 

2. for < t < T , for i — 1,2, for any 1 < r < +oo, 

\h{t)\r < |5o|re"- , \bi{t)\w^r < K\bo[u,,re~i 

3. forO<t< T, for i ^ 1,2, \ b,{t) < K' 

4. forO<t < T, fori = 1,2, 

mt)\l,2+Io\^^^\l,2 < eo 

mt)\Hs+J^\vn,\l. < K' 

5. forO<t< T, fori^ 1,2, 

\w^it)\l^^+J^\ym\l,^ < eo 
\m{t)\H^+J^\Vm\%. < K' 

then forO <t <T, 

\3w{t)\l,2 + C f \6w\l^, + \W{6w)\l^^ + \\^\{6w)\l,, 
Jo 

< c f ii + \wi\l,p + \yw^\l.) i\sb\l^p + \sh\l^^) .{61) 

Jo 

Proof. Combining (f62|l for i = 1, 2, we derive 

dtiSw) - (£ - a A) {dw) + ■ V) (Sw) - div (63 {V{Sw) + V^H)) 

= -(((5^7) . V)wi +div((J6)Vwi) (62) 
+ div(62V-Li?) + div(62V-L(55)) + div((5&)V-Lni) 

with n, R, Si and S2 obtained by solving 

div((l + 62)Vn) = div((l + 62)A((5w)-(i;2-V)((5i;)) 
div ({1 + 62)Vi?) = div f - a{v^ ■ V) (Sv) + (Sb) Avi ~ ((5^7) • V) wi) 
div((l + 6j)V5,) = div((l + &i)Awi - (i^i • V)wi) , fori = 1,2. 

Note that 02 = n + i? + ^2 and Hi = Si. 

Our strategy is again to muhiply (|62p by G^^{dw), integrate on and 
estimate each term arising to bound ^| Sw |^ 2- We deal with each term com- 
ing from the left member of equality (|62p as we did in the linearized vorticity 
equation (j4T|) . Let us only show how to deal with the other terms. 

First of all, let us emphasize that | C V(G''' /) ^ is controled by | V/ |^_2 + 
I ICI/ I', 2 provided that r + r' = -i. 
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First integrating by parts and applying Holder's inequalities, we obtain 
G-\Sw) div(S;i(5J7))| < \G^V{G-^{Sw))\2\G^iwi \p\61^\q 



/R2 

where 2 < q < +00 is such that ^ + ^ = 5 ■ Then using ([TT]) 

I / G-\Sw) div{w,iSP))\ < e\Giv{G-'iSw))\l + ^\w^\l^\Sh\l, 

(63) 

where e is intended to be chosen small enough. 
Similarly, with the same q, we also have 

G-\dw) div{{6b)Vwi)\ < \G^W{G-^{Sw))\2\G-^5b)\p\Vwi\q 



/R2 

thus using Sobolev's embeddings 

I / G-\Sw) div{{Sb)Vwi)\ <e\Gi V(G-i (Sw)) + - \ Vw, \5h\l^ 

(64) 

where e is again intended to be chosen small enough. 
In quite the same way, we obtain 

I / G-\5w)dW{b2V^R)\ < C|GW(G-i(,5w))|2|6oU,oo|Vi?|2 . 

Thus with pressure estimate (1^51) . estimate on Biot-Savart law (|17p . Holder's 
inequalities and Sobolev's embeddings we can derive 



1/ G-\Sw) div{b2\7^R)\ < e\Giv{G~'^ iSw))\ 



-\bo\l,oo\^w^\j,, \Sb\l 



G 



£ 

+ f l&o|?„,ool^i|^,2 (65) 
where e is once again intended to be chosen small enough. 

Using estimate p4| instead of estimate p3)) and inequalities ([T7| and p9|) , 

since ^ < i — i, we can obtain 
p I p' 

\[ G-\Sw) divib2V^ {5 S))\ < £\Giv{G-^{6w))\l 

JR2 

+ j\bo\loo\Sb\l 

x[\\/w,\j,,+\n,\l^^\w^\l^^]{66) 
where £ is still intended to be chosen small enough. 
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• At last, integrating by parts and applying Holder's inequalities, we have 
I / G-\Sw) div(((56)V^ni| < I V(G-i (Sw)) ^ | Sb V^Hi |, . 

where 2 < q < +00 is again such that ^ + ^ = ^- Again with pres- 
sure estimate (^5]) . estimates on Biot-Savart law pT)) and p^ . Holder's 
inequalities and Sobolev's embeddings we can derive 

1/ G-i((5w) div(((56)V^Hi| < £|G^V(G-i((5w))|i 




x[\\/w^\%,+\n\i^,\w^\l^p] (67) 
where e is once again intended to be chosen small enough. 
Gathering everything and integrating yield (|6ip . ■ 



4 Main results 

We now use our various estimates to derive our main results. 
4.1 Existence and uniqueness 

Before proving a result of existence and uniqueness of solution of equations ((TOl) , 
we state a lemma that will make a link between norm estimates and convergence 
of the iterative scheme. 

Lemma 3 Let T > and p > 1. Let (fk) be a sequence in L°°(0,r;R^), and 
(gk) a bounded sequence in L'^{0, T; R^) be such that for < t < T and k e N, 

fk+i{t) < fk9k ■ 
Jo 

Then (fk) is uniformally summable, namely, for < t < T , 

fc>0 

Proof. Using Holder inequalities and iterating yield fk{t) < K C'' , 
where G is a bound for (gj.) in LP{0, T; R+). ■ 

We can now prove the existence and uniqueness parts of Theorem 1. 
Proof. 



25 



• Existence. We build a sequence {{bk, Wfc))fceN* such that, for any k G N*, 
r 5t6fc+i + (K-ie)^V)6fe+i^ ^ = 

\ dtWk+i - (£ - a A) lUfe+i + (ufc • V) uj^+i = div [bk (Vwfe+i + V-'-IIa; 

where C and A are as in (jlip. (S^fc) is obtained from (wfe) by the Biot-Savart 
law, 

Vk — av""^ + Vk , Wk — aG + u'k , for k E N*, 
and (Vllfe) is obtained by solving, for any k E N*, 

div ((1 + bk) VUk+i) = div ((1 + bk) Avk+i - (vk ■ V) Vk+i) 

with initial data (bo^wo). For k — 0, we solve the system with Vk{t) = 
and bk{t) = 0. 

Let us show how we propagate bounds on {bk,Wk)- 

Step 1 Fix Kq > and choose eo > small enough. We can propagate 

1. for any 1 < p < +oo, 2 < q < +oo, thanks to Proposition 5, 

\bk{t)\p < \bo\pe-p , \bk{t)\^^g < |&oU,ge~ie^» 

2. and thanks to Prosposition 8, 

I Mt) \1.2 + Cko f (I Wk 1^,2 + I ^Wk 1^,2 + I lel^fe \l,2) 

Jo 

< CKo{\wo\l^2 + \bo\wA) 
which provides us, thanks to Proposition 1 and Sobolev's embeddings, 

\vk\8 < C\wk\s < C\wk\w,2 < Ko 

Io\^k\L < CJ^ ilwklC + l'^Wkll) < min(^,Xo) . 

Step 2 Again choosing Eq small enough independently of i,and using 
Proposition 2, we can obtain, when < s < 1 and 1 < s" < 2 — s, 

Jo\^Vk\m < C{t + tHlo\yvk\l.)'^) < Cit + j;;\wk\l.) 

\I'Vk{t)\2 < C\wk{t)\^,2 < CkoSo 

Io\I'vk\l.„ < C J^{\wk\l,, + \Vwk\l) < Ko 
and propagate, for < t < T, 

1. when 1 + s < s' < 2, thanks to Prosposition 6, 

I bk{t) 1^.. < Cko,T 

2. and thanks to Prosposition 9, 

\rwk{t)\l + C I \PVwk\l < Cko,t 
Jo 
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which provides us, for < i < T, 

' {\wk\l + \P^Wk\l)) < CK.n 



I iVwfcU.+i < C{t+ f { 
Jo Jo 



thus, thanks to Prosposition 6, for < i < T, 

\bk{t)\H^^2 < Cko,T ■ 

Step 3 Wc now prove the convergence of the scheme. Set {6b)k = hk+i — bk 
and {5w)k — Wk+i — "Wk- Choose niax(4, ^) < p < q. Propositions 7 & 
10 give us, for T > 0, for any < < < T, for any fc e N*, 

\mk+i{t)\i,, + \{sw)k+i{t)\i,, 

< CTj;,il + \wk\l,, + \Vwk\U {\{Sb)k\l, + \iSwh\l,2) 
for some < 77 < s such that ^ < p < +00. 

Now in order to apply Lemma 3 with fk = \ {Sb)k |^^p + | {6w)k \t,2> 
remark that 

- since (G-iwk) is bounded in L°°(R+; L2(r2)) and {V{G-iwk)) is 
bounded in (r+ ; ^2 ( ) ) ^ {G-^Wk )is bounded in L'' (R+ ; (R^ ) ) , 
for some 2 < r < +00, by interpolation and Sobolev embeddings 

- since (wk) is bounded in L°°{R+; L^{R^)) n L^iO,T; H''+'^(B?)), 
(Vwfe) is bounded in L'''(0, T; //''(R^)) for some 2 < r' < +00, by 
interpolation. 

Thus {bk) converges in L?^(R^) and (wk) in -^^(R-^), locally uniformally 
in t. 

This implies at once that (bk) and (wk) also converges in L^(R^). Now 
by interpolation 

— since (bk) is bounded in iJ^+^(R^), (bk) converges in (R^), for 
any < 77' < s + 2 

— since (wk) is bounded in 7J'*(R^), (wk) converges in (R-^), for 
any < 77" < s. 

These properties enable us to take the limit in the sequence of equations. 

Note that we recover the regularity on the limit by a mere application of 
Fatou's lemma. 



• Uniqueness. We obtain a bound similar to (|68[) for the difference of two 
solutions. Then Gronwall lemma gives the result. 
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4.2 Asymptotic behaviour 

We now state the asymtotic part of Theorem 1. Note that under the hypothe- 
ses of Theorem 1, the foUowmg asumptions are fulfilled. 

Theorem 2 Let a G R. For any < 7 < ^, there exist Eq > and K,K' > 
such that if {b,w) is a solution of ilO\) with initial data {bo,wo), such that 

I b() \w,2 < £0 , I ^0 U,oo < £0 J 

and for t>0, \ wit) + /J | Vw < eo , 

then, for t>0,\ wit) ^,2 < K' e'T* (| wq + | ^.,2). 
Proof. Let < 7 < 7' < i. Let us recall p5)) : 

j^G-^wCw< -7'|5;|2^^-(i_y)(i|V^5|^,2 + ^lf ■ (69) 

Then we deal with the other terms of the vorticity equation as we did to 
obtain estimate (j43p . except for the pressure term and 

I / G-'wv-Vw\<C\w U,2 (I w 1^,2 + I Vw 1^,2) (70) 

obtained thanks to estimate pS]) and Sobolev's embeddings. 
We treat the pressure term as follows 

1/ G^^w div(6V-Ln)| < C|G5 V(G-iw)|2|G"Hv-Ln|2 

with 

|G-Uv^n|2 < G(|&oU,2e-5|a(l + 6)Ai;G_a2(z;G.v)«G|^ 

+ I &o U,oo \il + b)Av-a ((w^ ■V)v + iv-V) w^) ~iv-V)v I2) 

and |(l + 6)Au|2 < G (1 + | 60 U,oo) | VSi I2 

|(i;G.V)?|2 < C\v''\oo\w\2 

\iv-V)v''\2 < ClWv^UlvU < C\w\^a 

\iv-V)v\2 < C\v\oo\yv\2 < Ci\w\^^2 + \Vw\2)\w\2 ■ 

This yields 

^^(1 ^ 1^,2) + 7l 55 \l,2 + Ci\Vw \l2 + 1^,2) <C\bo 1^2 e-* (71) 
which integrating gives our result, since 27 < 1. ■ 
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